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Abstract. We reconsider stability of the non-trivial fixed point in 6 — e dimensional 
effective action for the random field Ising model derived by Brezin and De Dominicis. 

C^ , After expansion parameters of physical observables are clarified, we find that the non- 

S ' trivial fixed point in 6 — e dimensions is stable, contrary to the argument by Brezin 

I ■ and De Dominicis. We also computed the exponents v and ri by the e expansion. The 

pH , results are consistent with the argument of the dimensional reduction at least in the 

Q ' leading order. 
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Tjj- ' 1. Introduction 

Q^ ■ Recently, Brezin and De Dominicis derived effective scalar field theory for the Ising 

Q^ ■ model in a Gaussian random field within the replica formalism ^. They showed that the 

^ • theory contains five cj)'^ coupling constants as well as the standard 0^ coupling constant. 

They also pointed out that we meet singular fluctuation on the critical surface below 
1-^ ■ dimension eight when we take the zero- replica limit. In order to resolve this problem, 

G . they redefined the coupling constants such that the beta functions for the new coupling 

O ' constants do not explicitly contain n, the number of the replica components. 

The effective theory contains not only the five 0^ coupling constants Mj, [i = 

S>^ . 1, . . . , 5), but also a parameter A with the mass dimension 2, which is related to strength 

of the random field. Since the new coupling constants introduced in reference |l| are 

defined as linear combinations of A-Uj, they have the dimension 2 + {4 — d) = 6 — d in 

(i-dimensional effective theory. 

According to the beta functions corresponding to the new coupling constants, there 
is a non-trivial fixed point in 6 — e dimensions as in the case of the pure 0"^ theory in 
4 — e dimensions. However, aspects of flow is different from the pure theory because 
some flow is repelled by the fixed point. Therefore, they concluded that the non-trivial 
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fixed point is unstable and the e-expansion cannot be carried out near six dimensions. 
It indicates that the dimensional reduction |^, Q breaks down near the upper critical 
dimension as well as the lower one @, 0]. 

However, it depends on choice of coupling constants in this model whether a fixed 
point is attractive or repulsive. Hence it is important to study how physical observables 
are expanded by Ui and A, and to know what the effective couphng constants suitable 
for perturbation are. To make this point clear, it is instructive to consider the following 
simple example. Let A be a dimensionless, renormalized coupling constant. Suppose 
that the beta function is given as 

(3x{\)=x\-yK\\ (1) 

where A = Am~^ and m is the renormalized mass, the inverse of the correlation length. 
We denote by x and y dimensionless coefficients. The fixed point A* is easily found as 
x/{yA). Then we have 

d(3x 



dX 



(2) 

A=A* 

Next let us introduce a new coupling constant by 



g ^ AA. (3) 

The beta function for g is computed to be 

Using (H) and (§), one finds 

^ g=9* 
where g* is the zero of Pg. Suppose that < x < 2. Equation (Q) shows that the fixed 
point is unstable. On the contrary, it seems stable from the result @. This example 
shows that we have to clarify which coupling constants play a role of an expansion 
parameter for studying the stability. Note that the existence of the second term in the 
right-hand side of equation (^) shows breakdown of the covariance of the beta functionQ. 
It appears because the fixed dimensional parameter A enters the beta function and the 
definition (^. 

Although Brezin and De Dominicis chose Aui for all i as the coupling constants 
in order to eliminate n dependence, it has not yet confirmed whether the new coupling 
constants correspond to good expansion parameters. The main purpose of this paper 
is to obtain the good expansion parameters and to reconsider the instability studied 
in reference [Q. We conclude that their beta functions are not for the good expansion 
parameters. 

This paper is organized as follows: in the next section we perform the naive 
perturbative expansion of the free-energy density in the disorder phase and determine 
the expansion parameters. Thanks to a finite mass in the disorder phase, there are 
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no singular behavior in the hmit of n — ^ discussed in reference |I|. Therefore we 
can take the zero-rephca hmit at the very beginning. We find that the good expansion 
parameters are obtained from the effective couphng constants A^'-Uj, where ai depends 
on i (see equation (0)). Our effective couphng constants have various mass dimensions 
at the tree leveL One of them becomes marginal but the others are irrelevant in six 
dimensions. When the dimensions lower by e, the marginal coupling constant becomes 
irrelevant at the non-trivial fixed point due to the one-loop correction, while the other 
coupling constants remain irrelevant. 

In section 3, we turn to renormalized perturbation, where we define the 
renormalized effective coupling constants and compute their beta functions. Computing 
the scaling matrix (i.e., the derivative of beta functions at a fixed point), we find that 
the fixed point is stable. 

The correlation-length exponent and the anomalous dimension of the correlation 
function in c? = 6 — e are computed in section 4. The results are identical with the case 
of the d = 4 — e pure 0^ theory in the leading order of the e-expansion. It is consistent 
with (d, d — 2) correspondence 0, |^. 

The last section is devoted to summary and future problems. 

2. Coupling dependence of the free energy 

Here we observe the bare-parameter dependence of the free-energy density. 

The (i-dimensional effective action derived by Brezin and De Dominicis is given by 

O = Do + Ointj 

a, 13=1 
^int = J d X i^ — a^ + -^(^3^1 + Y^2 + ^^2^1 + -^cT^j , (6) 

where a and /3, which run from 1 to n, specify a replica component and 

n 

a,^Y.iMx))\ (A: = l,...,4). (7) 

a=l 

The interaction conjugate to Ui contains one summation over the replica index while the 
other interactions have more than one, which apparently shows that the other coupling 
constants M2, • • • , M5 are less relevant in the zero- replica limit. 
The quadratic part Sq defines the propagator: 

^°/5(^) = ^T^*^"^ ^ {p' + ty + °^''^' 

= A{p)6^p + Bip) + 0{n). (8) 

Note that B(p) dominates over A{p) in low-momentum region with small t, while A(p) 
can be more relevant in the replica limit because 6ai3 associated with it reduces powers 
of n. 
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A simple dimensional analysis gives that the parameters of this theory have the 
following mass dimensions 

H=4-d(^ = l,...,5), [t] = [A]=2. (9) 

We note that equation (H) indicates that the correlation length ^ is proportional to t^^/^ 
at the tree level. We can see the critical phenomena by one-parameter tuning t = 
despite that A is a relevant coupling constant. 

Let us look at the perturbative expansion of the free-energy density /, which is 
obtained as 

/ = lim :^, e-^^-- = / 17 I^^e"^, (10) 

n^O n J \ 

where V is the volume of the system. The leading order is easily computed. The result 
is written as 

Here A is the ultraviolet cutoff. The integration is identical with the free-energy of the 
Gaussian model, whose singular part behaves as f^^^ when t is sufficiently small. Hence 
the singular part of / is 

fi^^o,M'"^-\l + 0{^-H)), (12) 

where the O(A^^t) term is irrelevant near the criticality: t ~ with fixed A. 
Next, we compute the first-order correction 

limi-(5i„t). (13) 

For instance, the Ui interaction gives 

^(^i + A^2)', (14) 



A= /tt^T^^TI' ^ = 1'2- (15) 



where 

d'^p 1 
(2vr)^ {p^ + 1) 
Extracting the singular part of Ak from equation (p!^ , we have 

!^ (a,trf/2-i + Aaat'^/'-')' = ^A^a^t'^/^-^ (i + Q (A-^t)) , (16) 

where Oi and 02 are some constants. Similarly, collecting the leading contribution from 
the each vertex that survives in the zero- replica limit, we get 

/,ing = At 2-1 lao+{ ^alt-^~^ + '^ ' aia2t"^ + -—Aalt'^-^ \ \ , (17) 

up to the first order of Ui. Similar computation to the first few loops implies that the 
expansion parameters are A"'Ujt'^/^~^~"\ Here 

ai = 1, 0^2 = ^3 = 0, 04 = —1, as = —2. (18) 
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We can see that the Aui term survives while the other corrections in (|T^ vanish 
in 6 — e dimensions near the criticahty. 

Do our expansion parameters work for all diagrams in /smg? Let us consider an 
arbitrary connected diagram in /. It contains two kinds of internal lines, A{p) and B{p). 
Let ^A and ^B be the number of A{p) and B{p) in the diagram respectively. Since 
B[p) is proportional to A, the diagram generates the following correction: 

F(t,A)A#^n«f, (19) 

1=1 

where F(t, A) represents momentum integrations and summations over the replica 

indices. We find from the dimensional analysis that ( |T9t) has the following form: 

F(l, A/t) Af^/^-i I^^-HY' n (A^'w.t^/'-'^"')"^ . (20) 

i=l 

Comparing ([19D and (^) we find 

Po= \y^<^^P^+^-*B. (21) 

We must show that 

po > 0. (22) 

If Po < 0, the expansion may fail because the factor At~^ grows near the criticahty, even 
though the last factor gives small correction. 

Here we prove (l2^). Because the number of the internal lines are 2Y^pi, 

#^ + #5 = 2Ep- (23) 

i=l 

From equation (|]) we find that the interaction conjugate to Ui has 2 — Oj summations 
of the replica indices. Thus the number of the sums in F is 

5 

E(2-«.)P^- (24) 

It should be reduced up to 1 by Kronecker's delta associated with A{p): otherwise, F 
becomes O(n^) and disappears in equation (|l^. Thus, we have the following restriction 

to#A: 

5 

#A>E(2-a,)p,-l. (25) 

i=l 

Using (|23| ) and (pSf) , we obtain (|^). 
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Figure 1. Feynman diagram for the beta functions (E9). 

3. Renormalized perturbation 

Now we are going to the renormalized perturbation. Using the renormahzed 2-point 
vertex function T2af3{p)i we introduce the renormahzed parameters m and Ar as 

r2 ap {p) = if + m^) 5^p + Anm^ + O (p^) . (26) 

We have shown that the expansion parameters are A^'Ujf^'^^^""', so that we put the 
following renormalization prescription by using the renormalization counterparts for 
them: 

A°T4,(0K-^ = ^„ (^ = 1,...,5), (27) 

where r4j(p) is the renormalized 4-point functions associated with u-i. Then one can 
perform the renormalized perturbation with the small parameters gi. 

The one-loop beta functions for gfj in 6 — e dimensions are easily computed by the 
usual manner. Following reference |I[], we keep the terms proportional to 

d'^p Ar 
(27r)^(p2 + m2)3' ^^^'^ 

which are expected to be important in 6 — e dimensions. See figure |l|. Here an internal 
line with a cross(x) represents B{p) and without crosses A{p). 
Similar calculation performed in reference [|l| gives 

/3i= -tgr + ^Kgl 
p2 = {2-e)g2 + 3Kgi{g2 + g3), 
p3 = {2-e)gs + 5Kgi{g2 + g3), 
(3, = iA-e)g, + AK{g2 + g3)^ 

/35 = (6 - e)g, + mK{g2 + (73)^4, (29) 

where K = r(4 — d/2)/{ATxY^'^. The non-trivial fixed point is, up to 0(e), 

^i* = ^e, gl = ...gl = Q. (30) 

Now we compute eigenvalues of the scaling matrix 



% 



(31) 
9=a* 
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The result is 

{e, 2-e, ^(3 + e),4-e, 6 - e}. (32) 

The eigenvalues are all positive. Thus the non-trivial fixed point is stable, contrary to 
the conclusion in reference [Q. 

The discrepancy stems from the ways of redefinition of the coupling constants. 
Brezin and de De Dominicis have essentially chose Auj as the effective coupling 
constants. The linear terms of the beta functions correspond to Auj are — e for all i, 
and negative eigenvalues emerge. However, the effective coupling constants that appear 
in perturbative expansion are A"*Uj as we have seen in the previous section. Since A 
has the finite mass dimension, scaling behavior of Auj does not refiect the stability of 
the fixed point. 

4. The e-expansion near six dimensions 

Finally, we compute the leading order of the e-expansion for the critical exponents v 
and rj at the non-trivial fixed point (|5(]|). 

To compute the critical exponents one can take account of the leading infrared 
divergence only. Figure |^ represents the most singular Feynman diagrams contributing 
to the self-energy. One easily sees that only the Ui term contributes to the most 
singular diagrams; vertices M2,...,M5 combined with B{p) giving additional n factor 
are suppressed by the ra — > limit. Thus, from the renormalization conditions (^) 
and (pTl) , the one- loop correction of the bare mass and the two- loop computation of the 
field renormalization factor Zi are written in terms of the renormalized parameters as 
follows: 

From equations (^91) and (|33D , the composite field renormalization factor Z2 is written 
at the one-loop level as follows: 

d^'^ . ir(3-t//2) ^, 2. 



dt 





Figure 2. Feynman diagrams contributing to the self-energy. 
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Following the standard argument B H], define 



72(5'!, •••,fl'5) = m— InZa 



A,Wi:fixed 

(36) 



AjUiifixed 

We can compute the critical exponents v and rj from the following relationship: 

1 



V 



V = li{9*i,---,9*5)- (37) 

Substituting the values of the fixed point (|30|) into the above equations, we obtain the 
leading order of the e-expansion for the critical exponents u and r], 

These are exactly the same as the leading-order results in the d = 4 — e pure (p^ theory 



5. Summary and discussion 

In this paper, we have clarified effective coupling constants for perturbative expansion of 
physical observables in effective theory for the random field Ising model and reconsidered 
the stability discussed in reference [0. Our result ensures the e expansion near the 
upper critical dimension. We have also calculated the exponents z/ and 77. The 
result is consistent with the Parisi-Sourlas dimensional reduction in the leading-order 
computation. It is a nontrivial problem whether the consistency is preserved beyond 
the leading order, which will be reported elsewhere. 

Our beta functions (^) tell us that the upper critical dimension is six, which is 
consistent with the rigorous result [§. Further, they show that g2 and 93 change to 
relevant coupling constants when d becomes lower than four. Thus it does not surprise 
us that the dimensional reduction does not hold when d = 3 M- 
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of the manuscript . 
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